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O ! Abstract 

J>^' We review the theory of hght propagation in moving media with 

1-5 ' extremely low group velocity. We intend to clarify the most elemen- 

tary features of monochromatic slow light in a moving medium and, 
whenever possible, to give an instructive simplified picture. 



1 Introduction 

Waves experience moving media as effective gravitational fields. To under- 
stand why, imagine a definite example — light traveling in a moving trans- 
parent fluid such as flowing water |l|, §]. In each drop of the liquid, light 
propagates along a straight line and so each drop distinguishes a particular 
inertial frame. Now imagine light passing from one drop to a next one which 
happens to move with a different velocity vector. Again, the new drop in the 
way of the light distinguishes an inertial frame, but this new frame will differ 
from the previous one. Light traveling in a non-uniformly moving medium 
is constantly forced to adapt to new inertial frames. An analogous situation 
occurs in curved space-time when light, and all other matter, experiences 
space and time as consisting of connected inertial frames. Space-time can be 
thought as being locally fiat, even in the vicinity of the most violent gravitat- 
ing objects, yet each local frame is non-trivially connected to the neighboring 
frames, a connection mathematically described in terms of the curvature ten- 
sor 0. In a moving medium, the local co- moving frames of the medium play 
the role of the local pieces of fiat space-time. Because theses pieces differ 
in non-uniform flows, light propagation in a moving medium resembles light 
propagation in curved space-time, i.e. the medium appears to light as an 
effective gravitational field @l-0. 

One could conceive of employing moving media for creating artificial as- 
tronomical objects in the laboratory. For example, water going down the 
drain of a bathtub appears to light as a rotating black hole. However, in 
order to observe spectacular effects of general relativity in laboratory-based 
analogues, truly astronomical flow velocities are required. For establishing 
a black hole, the fluid should move faster than the speed of light in the 
medium (c divided by the refractive index n). The chances of creating anala- 
goes of black holes are much better when one employs sound instead of light 



in appropriate supersonic flows 0-|16]. One of the most fascinating effects 



of black holes is Hawking radiation [|r^-[|iy], the spontaneous generation of 
photon pairs near the hole's horizon where one of the photons falls into the 
hole and the other tunnels out of the attraction zone and become visible. The 
acoustical equivalent of Hawking radiation is quantum sound in superfluids 



[0]-|]T6|. However, superfluidity tends to break down before the fluid has a 
chance to move faster than the speed of sound [^ . Furthermore, quantum 
sound and flow are just two aspects of the same object, the quantum liquid. 



The flow is the macroscopic and the sound the microscopic motion. Under 
extreme circumstances such as near the horizon of a sonic black hole, a clear 
distinction between sound and flow might be difficult. An advantage of light 
in a moving medium is the clear separation between wave and flow, and both 
light and medium can be regarded as separate quantum systems. 

Recently, light has been slowed down dramatically [^-p3[ due to an ef- 



fect called Electromagnetically Induced Transparency [^]-[26|. The use of 
slow light could open the opportunity for observing the radiation of quantum- 
optical black holes. Note, however, that slow light is a more complicated 
phenomenon than light in a moving non-dispersive dielectric where the re- 
fractive index does not depend on the frequency of light. Slow light is based 
on a highly dispersive medium created by dressing the atoms of the medium 
with an appropriate light beam. Instead of the phase velocity only the group 
velocity is reduced. Therefore, it is not always advisable to conclude directly 
from the behavior of light in ordinary dielectrics to the motional effects of 
slow-light media. It is the purpose of these notes to clarify the most elemen- 
tary features of slow light in a moving medium [^ and, whenever possible, 
to give an instructive simplified picture. 

2 Dispersion relation 

Imagine the medium as being decomposed into drops. Each drop should be 
small enough such that the flow does not vary significantly within the size of 
the drop, but each drop should be sufficiently large to sustain several optical 
oscillations. We thus assume that the wave length of light is much shorter 
than the typical scale of changes in the flow. In this case, we can describe 
light by a dispersion relation in the local co-moving frame of each medium 
drop denoted with primes 

*'^-^-x(^')i^ = o. (1) 

The susceptibility x{^') characterizes the medium. In Electromagnetically 
Induced Transparency |2^-|2^ a coupling beam dresses the upper two levels 
of a three-level atom shown schematically in Fig. 1. The dressing of the 
upper two levels influences strongly the propagation of the probe light with 
a frequency that should match the atomic transition frequency ujq between 
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Figure 1: Three level system needed for the creatfon of electromagnetically 
induced transparency. A strong drive laser field couples the levels a and 
c, making the medium transparent for the weak probe laser tuned to the 
transition a^^b. 



the lower and one of the upper levels. Exactly on resonance, the medium 
decouples from the probe light and becomes transparent, whereas ordinary 
dielectrics are extremely absorbing at atom-light resonances. Furthermore, 
the susceptibility changes rapidly, see Fig. 2, and, near the resonance, x{^') 
assumes a linear dependence on the detuning between the light frequency u' 
and the atomic frequency uq, 

Xiu;') = ^{uj'- uo) , (2) 

UJoVg 

giving rise to a very low group velocity Vg. If the slow-light medium is moving 
coherently, the motion will slightly detune each atom from exact resonance 
due to the Doppler effect. 

Imagine that the flow of the medium, u, varies only in two dimensions 
and that the coupling beam propagates orthogonally to the flow, see Fig 
3. In this case we can ignore the flrst-order Doppler effect of the coupling 
beam and we should focus only on the Doppler detuning of the probe. In the 
laboratory frame the probe shall be monochromatic at the atomic frequency 
ujq. In the co-moving local frames of the medium the probe-light frequency 
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Figure 2: Susceptibility for the probe laser beam in electromagnetically in- 
duced transparency. The plot shows the dependence of the real (whole lines) 
and imaginary (dashed lines) parts of the susceptibility x on the detuning A 
of the probe beam. The drive beam is assumed to be on resonance. Arbitrary 
units are used. 



is Doppler-shifted with the detuning 
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to first order in u/c. Let us transform the dispersion relation (|l]) to the 
laboratory frame. We take advantage of the Lorentz invariance of k"^ — {oj/cY 
and obtain to first order in u/c the relation 
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with 



ko = — 
c 



(5) 



In Ref. [^ a more complicated dispersion relation has been derived which is 
correct up to second order in u/c. Note that the relation (§) describes effects 
in leading order of u/{vgcY^'^, because the second-order corrections of Ref. 



27] are proportional to v~ 




Figure 3: Slow light in moving media. The medium should move orthogonally 
to the coupling beam (grey arrow). This beam prepares the medium such 
that the group velocity of light is significantly reduced. The probe beam 
(lines) experiences the Doppler detuning due to the moving medium and 
becomes deflected or even trapped. 



3 Magnetic model and metric 

We are going to use the simplified relation (^ to illuminate the most char- 
acteristic features of slow light in moving media. First we note that we can 
write the relation as 

{k + k,^J-kl^^=kl. (6) 

Monochromatic light waves with fixed polarization thus obey the wave equa- 
tion 

L^^ + ko-] <P-kl-^<p = kl<p. (7) 

The flow has a two-fold effect: On one hand, the velocity u appears as an 
effective vector potential, for example the magnetic vector potential acting 
on an electron wave [^], and, on the other hand, the hydrodynamic pressure 



proportional to u ||2^ acts as a scalar potential. This resembles the Rontgen 



effect of static electromagnetic fields on polarizable atoms [^0|, pl |. 

Alternatively, we can regard the flow as generating an effective gravita- 
tional field. To see this, we introduce the four-dimensional wave vector 

k, = (/Co, -k) , (8) 

and, adopting Einstein's summation convention, write the dispersion relation 
©as 

g^'^k.k, = (9) 

with the contravariant metric tensor 

«- - ( u/\, "-'? ) • (1°) 

Four-dimensional wave vectors are null vectors with respect to the effective 
metric g^^ and, in turn [||], light rays follow zero-geodesies with respect to 
the line element 

ds^ = g^ydx^dx'^ , dx'^ = {ct, x) . (11) 



The covariant metric tensor g^i, is the inverse of the contravariant one, and 
is given by 



In contrast to sound |jTj] or hght in non-dispersive media |]l[, a genuine 
event horizon of slow hght cannot exist at this level of approximation. The 
potential existence of an event horizon is thus confined to effects to second 



order in M/(fgc)^/^which are excluded from our simplified model [^ 

4 Non— relativist ic analogue 

To get more insight into the behavior of slow light in moving media, let us 
study another analogue. Consider a non-relativistic free particle attached to 
moving frames. The particle is supposed to move freely in each of the local 
frames but is bound to adapt from one frame to the next, similar to light 
in drops of flowing water. However we use pure non-relativistic physics to 
transform from frame to frame. In case the local frames constitute a global 
frame, for example a rotating solid body, the particle would move along a 
straight line with respect to this frame. In the laboratory frame, the parti- 
cle's trajectory appears to be bent, attributed to the effect of Coriolis and 
centrifugal forces. The Lagrange function of such a fictitious non-relativistic 
particle is 

L = \v" = \{^-^)\ (13) 

up to a constant mass factor that we can set to unity, because we are inter- 
ested in inertial effects only. The velocity of the local intertial frames is w 
and we have used the non-relativistic addition theorem of velocities 

v' = V - w . (14) 

Note that we could apply Einstein's relativistic addition theorem of velocities 
to describe light in moving non-dispersive media of refractive index n, but we 
should use an effective speed of light of c/n in the theorem. The momentum 
of the fictitious non-relativistic particle is 

dL , , 

p = — = V - w (15) 



and the Hamiltonian is 



H = v-p-L = -p'^-w-p = E. (16) 



The important point is that we can translate the Hamiltonian ([T6| ) into the 
dispersion relation (H) by setting 

p = c— , w=-u, E = -c\ (17) 

ko Vg 2 

with the enhanced effective frame velocity w [^]. Consequently, slow light 
experiences the moving medium in the same way as a non-relativist ic particle 
experiences moving inertial frames. Excused by its fictitious nature, the 
"non-relativistic" particle would move at an extraordinary speed reaching c 
in regions where the medium is at rest, because the energy is 

(18) 



(19) 

UL 

the equation of motion of the fictitious particle 

v = (Vxw)xv+- Vw^ . (20) 

The first term describes the Coriolis and the second the centrifugal force. 



5 Light rays around a vortex 

Consider a specific example of a non-uniform flow, a vortex. Quantum vor- 
tices in alkali Bose-Einstein condensates have been recently made |^, P5| ]. 
They give rise to intriguing slow-light phenomena. However, vortices do not 
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obtain from Hamilton's equations. 
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dp 
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posses genuine event horizons |32| and they will not radiate spontaneously. 



The flow of an ideal vortex of vorticity W is 



u=—e^ = i—, (21) 



written in polar coordinates or using the complex notation z = x + iy, z* = 
X — iy oi the planar Cartesian coordinates x, y. The curl of the vortex flow 
vanishes, 

V X u = , (22) 

and hence we obtain the equation of motion of a particle in an attractive r~^ 
potential |]36|-|]38 

,2 



^V..'^-^i, (23) 



9 



or, in complex notation. 



^ + —^Tzr.=^- (24) 



t;? \z\^ 



The trajectories of the fall into a r ^ singularity are well known in polar 
coordinates |^ but here we suggest a more elegant description using the 



complex notation of the planar Cartesian coordinates, 

.= (^-c« + ,-j [^^ct~,-j (25) 

with the real parameter b and 

/ y^2 \ -1/2 



One can easily verify that this solution satisfies the equation of motion (|2 
In the infinite past, t — > — cx), the trajectories approach 

z ct + ib . (27) 

From this asymptotics we see that the light is incident from the right to the 
left with initial velocity c and impact parameter b. Depending on the impact 
parameter, the coefficient fi is real or purely imaginary. In the case of a real 
fj, when W^ < Vgb"^ the incident light particle will be able to escape, because 
the modulus squared of z, cH'^ + b"^ / fj,'^ , will never reach zero, the vortex core. 
In the infinite future, t -^ +oo, the particle will approach 

z r^ ct {-lY + ib = cte'''^' + ib (28) 

10 




Figure 4: Light rays spiraling towards the vortex core. The plot shows the 
hmiting spiral of described in Eq. (|2l 



and, consequently, leaves at the angle vr/x. In the case of a purely imaginary 
coefficient fi when W^ > f|6^ the light is sucked into the vortex, because the 
modulus squared of z, c^t^ + 6^//i^, reaches zero at the time t = — |6/(/ic) |. In 
practice, the light will hit the vortex core and will bounce back. In theory, the 
fictitious light particle disappears from this world and reappears on another 
Riemann sheet. In the borderline case when W^ = Vgb'^ the parameter /^ 
tends to infinity and the trajectory approaches 



, f 1 ibY 

z = —ct lim 1 

M^°o \ 11 ct I 




exp 



'VgCt^ 



(29) 



Light falling into the vortex core describes distinctive spirals illustrated in 
Fig. 4. 

Light waves around a vortex 

Near the vortex core, rays of slow light feel the presence of the centrifugal 
r~^ attraction but no Coriolis force, because the curl of the fiow vanishes. 
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Light waves, however, will show a distinctive interference pattern, an optical 



Aharonov-Bohm effect |^-|^ due to the long-range nature of a vortex flow. 
Far outside the vortex core, light rays are hardly deflected, yet rays passing 
the vortex in flow direction will be dragged and those swimming against the 
current lag behind. This does not affect the ray trajectories but it produces 
a phase shift between the dragged and the lagging rays, and, in turn, creates 
a typical interference pattern |41, 42 1. Note that the Aharonov-Bohm effect 



of waves in moving media is not restricted to light. Indeed, Berry et al. 



4^ have reported and analyzed a beautiful bathtub experiment with water 



waves in a tank, visualizing clearly the interference structure generated by 
the Aharonov-Bohm effect. A further experiment has been reported [^ 



showing spiral waves. Acoustical analogues of the effect have been observed 
in moving classical media ^^ and are predicted for superfiuids |^. The 
acoustical effect might even lead to a friction felt by traveling vortices due to 
the so-called lordanskii force [0-[^. The optical Aharonov-Bohm effect of 
slow light [^ can be applied to observe in situ the flow of quantum vortices 



in alkali Bose-Einstein condensates p^ , p5|] using phase-contrast imaging 



Consider the propagation of slow-light waves through a vortex flow. In 
polar coordinates the wave equation (^ is 



92 



1 d 1 



/'jiY' ^ iy* /'j'Y' 



'd_ 

dip 



+ iVAB + 



^AB 



kl 



(30) 



with the optical Aharonov-Bohm flux quantum 



W 



T^AB 



(31) 



We could simply eliminate the optical analogue of the vector potential by 
representing as 

= 00 exp(-zz/ABV') (32) 

where 0o feels only the local r~^ attraction due to the centrifugal force. The 
modulation exp(— iz/A^i/)) describes the long-range Aharonov-Bohm phase 
pattern |^T|. Strictly speaking, however, the ansatz (^) is only justified when 
the function exp(— iz/^BV^) is single-valued after a complete cycle of (^, i.e. 
when vab is an integer. For non-integer flux quanta the interference pattern 
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is slightly more complicated, showing, most prominently, a line of zeros of the 
wave function after passing the vortex (which resolves the problem of a multi- 
valued phase, because a phase is not defined when the amphtude is zero). 
Even in the case of non-integer optical flux quanta i'ab, the modulation 
expi^—iuAB^) accounts for the dominant phase pattern. 

The correct positive-frequency component of a slow-light wave incident 



from the right is an appropriate superposition |^| 



E 



m=—oo 



of the partial waves 

4>m = Um{r) exp{irrnp - iuot) , Um = {-i)" Ju{kor) 
characterized by the Bessel functions Jy with the index 



V 



m + VAsf - ^. 



AB 



(33) 



(34) 



(35) 



Figure 5 shows plots of slow-light waves around a vortex to illustrate the 
long-range phase shift brought about by the optical Aharonov-Bohm effect. 

Close to the vortex core, a light ray is doomed to fall into the singularity 
when the square of the impact parameters, 6^, does not exceed W^/f^. In 
this case the angular momentum m lies in the interval 

- 2vAB < m < , (36) 

because, according to Eqs. ([T5|jr7|) and the pattern (^) of the vortex flow, 

r X k = r X A;o(v - w)/c = {k^b - vab) ^z ■ (37) 

The fall of light rays into the vortex core corresponds to a drastic change 



of the behavior of the corresponding partial waves </)„. In the case (|36|) the 



index (^5[) of the Bessel function in a radial wave ( |3^ is purely imaginary. 
Using the first term in the power-series expansion of the Bessel functions 



nP we see that Um is rapidly oscillating near the core 
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Ur 



T(u+1] 
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exp 



ilm z/ In 




(38) 
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Figure 5: Slow-light waves traveling through a vortex. The dragging effect 
of the vortex shifts the phase of the incident light, depending on whether the 
light propagates with or against the flow. 
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This describes a rapid flow of liglit towards tlie centre of attraction, if Im u 
is negative and an outgoing wave for the positive branch. Far outside the 



vortex core, r 



Ur. 



oo, the radial waves approach 



V^vrAJor 



exp ( ik^r 



mu 



«^ ) +exp 



-ik^r + i 



n 



(39) 



as we see from the asymptotics of the Bessel functions |51|. For a negative 
imaginary u the second term in Eq. ( pQ]) dominates and so most of the light 
is incident, yet a small fraction with weight exp(— 2i7ri/) is able to tunnel out 
of the attraction zone. 

Where is the borderline between pure influx and partial tunneling out of 
the attraction zone? Consider the geometrical optics of a radial wave. We 
make the eikonal ansatz 



Ur 



'^m^m ~r '^m^m 



u- 



\um\ exp(±z_R) 



(40) 



regard \urn\ as a slowly varying envelope, and obtain from the wave equation 
30f) the radial eikonal 



R 



IV 



\ p^ + 1 — arsinh - + i— 



P 



iv 



(41) 



Quite typically, valuable insight into the behavior of waves can be gained 
by analytical continuation to the complex plane of the wave's argument. 
Consider the line in the complex plane of the radius r where the eikonal is 
purely imaginary, i.e. where 



Re 



p2 + 1 — arsinh — 
\P. 







(42) 



see Fig. 6. This line is called a Stokes line |£2|, Q. When passing a Stokes 
line, one component of the superposition (|40| ) becomes exponentially small 
and hence irrelevant |5^ . The asymptotics in the limit of geometrical optics 
switches from a two-component regime, for example in-and-out flux, to a 
single component behavior, such as pure influx |]S2|. The Stokes line in Fig. 
6 thus marks the boundary where light is able to tunnel out of the optical 
black hole. 
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Figure 6: Stokes lines of an optical black hole. Stokes lines |^ are lines in 
the complex plane of the radius r where the radial eikonal of light waves is 
purely imaginary. The lines connect the turning points of light rays which 
are complex in the general case of trapping (and purely imaginary for an 
optical black hole). The Stokes lines are calculated from the radial part of 
the eikonal (41) and are displayed in the scaled units employed there. The 
arrows indicate the light flux in the regions bounded by a Stokes line. 
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6 Summary 

We analyzed the propagation of slow light in moving media in the case when 
the light is monochromatic in the laboratory frame. Slow light is generated 
by the electromagnetically induced transparency of an atomic transition. 
The Doppler detuning due to the moving medium generates flow-dependent 
effects. We considered analogies to magnetism and general relativity and 
studied light propagation around a vortex in some detail. The subject is far 
from being exhausted — light detuned from the atomic resonance and non- 
monochromatic light will behave differently, opening perhaps new exciting 
analogies and opportunities for laboratory tests of astronomical effects on 
Earth. 
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